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Einleitung 

Tullio Levi-Civita hat die fundamentale Erkenntnis gefunden, daß die 
Gleichungen für die Bewegung der Planeten auf die homogene Schwin¬ 
gungsgleichung zurückgeführt werden können. Diese Erkenntnis wurde 
von Pythagoras und Kepler schon in der Form erahnt, daß den Planeten 
Töne zugeordnet wurden. 

E. Stiefel hat das Werk in seinem mit G. Scheifele verfaßten Buch [7] 
dargestellt und weitergeführt. Die Erkenntnis von Levi-Civita ermö¬ 
glichtes nun, dieTheorie der Gleichverteilung auf dieTheorie der Plane¬ 
tenbewegungen anzuwenden. 

Es ist weiters bekannt, daß auch das Modell von Friedmann-Lemaitre- 
Milne-Heckmann, und damit die Friedmannsche Gleichung auf die 
homogene Schwingungsgleichung zurückgeführt werden kann. Die Tat¬ 
sache, daß es mehrere Planetenbahnen gibt, hat mich dazu angeregt, rein 
mathematisch auch mehrere Welten einzuführen. Damit will ich nichts 
über die tatsächliche Existenz solcher Welten aussagen (Dies gilt viel¬ 
leicht noch mehr für die Anmerkung 4). Es ist mir bewußt, daß ich mich 
damit am Rande des für einen Mathematiker Erlaubten bewege. 

Als Ergänzung zu der im Text angegebenen Literatur verweise ich auf 
das Werk von Jan von Plato [6]. 
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Ich möchte noch bemerken, daß die Kenntnis von Teil I nicht voraus¬ 
gesetzt wird, die Lektüre des Epiloges wird allerdings empfohlen. 

Zum Schluß möchte ich Frau Marianne Wenger meinen herzlichsten 
Dank für die Abfassung des Manuskripts aussprechen. 


1 


Wir betrachten zunächst die Gleichung 


d 2 r K 


(i) 


nach der Methode von Levi-Civita, die dann vor allem von E. Stiefel ([7]) 
entwickelt wurde. 

Wir gehen direkt vor, indem wir die von Levi-Civita angegebene 
Lösung verifizieren. 

Es sei tu eine beliebige positive Zahl, und wenn K > 0 ist, dann defi¬ 
nieren wir die positive Zahl a durch 


a 


2 


K 

2 üj 2 


( 2 ) 


Dann setzen wir als Lösung von (1) an: 

r = a 2 (cos(ü/r + 'ip)) 2 (3) 


wo ip modulo 27t bestimmt ist und definieren die Zeit T durch das 
Integral 


/ = 



( 4 ) 


welches nur bis auf eine additive Konstante bestimmt ist. 
Es ist also 

dt 

dr 


Wir nennen r die fiktive Zeit im Gegensatz zur Zeit t. 

Wir zeigen nun, daß die Funktionen r(r) und t(r) eine Parameterdar¬ 
stellung der Lösungen von (1) sind. Zu diesem Zweck definieren wir als 
Geschwindigkeit 



( 5 ) 
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und benützen, daß 


ist, und erhalten 


dr I dt 
dr / dr 


sin 2 (üJT-\- , ip) , 

v = 2 cj T~, - 7ö7 -Wm = utgaj(r ff ip). 

1 ff cos(2 (COT + ip )) 

Es wird also (wir setzen rj = tur ff ip) 


V 2 K . J, , 2 1 


— -~ = 7w (tg r/j - 


2 r 


(cos rj) 


-u>‘ 


( 6 ) 


da 


ist. 


— y~ ~ ’n = 1 

COS z Tj 


Es ist nun 

dv d 2 r 
dt dt 2 ’ 

die Differentiation von (6) ergibt (1). Dabei muß 

dt 


dr 


/0 


sein. Es ist rückblickend (6) die Energiegleichung mit der Integrations¬ 
konstanten uj 2 Es ist 


u = a cos(u/r ff t/>) 

die Lösung der Schwingungsgleichung 

d 2 u 


dr‘ 


ff tu u — 0 


( 7 ) 

( 8 ) 


mit den Anfangsbedingungen 

#(0) = a cost/i, u ) = — cjrzsin , 0. 

Wählen wir nun eine gleichverteilte Folge ((/ff) modulo 1, so erhalten wir 
die Schar der Lösungen (t/ff = 2i r</ff) 

Uk(r) = a cos (tur ff t/ff), 


( 9 ) 
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und für (1) die Schar der Lösungen 

r k (r) = u\{t) = ü 2 cos 2 (urr + ip k ) 

und 


tk I f'k(o' s jdo' H - Cki 
Jo 

wo Ck Konstante sind. 


( 10 ) 


Wir betrachten nun s solche Gleichungen (l 7 ) mit Konstanten K J , also 

d 2 E K J 


+ 


dt 2 (rJ) 


0 


0 ') 


(j = 1,... ,s). Ein Beispiel: Es seienr Teilchen p J mit den Massen m J für 
j = 1,..., s gegeben, die von einem Stern S mit der Masse M nach dem 
Newtonschen Gesetz angezogen werden. Dabei ist dann 

K J = G(M + m J ) 

(G > 0 Konstante). 

Es seien die zugehörigen Energiekonstanten tu 2 ,. , c uf mit 


die Gesamtenergie. 
Weiter sei (aj positiv) 


üü 2 + + u) 2 — E 


KJ 


J 2tü 2 ’ 


( 2 ') 


dann sind wieder Lösungen des Systems (l 7 ) gegeben durch ('ipj =2 tt ipj) 

r k = U jk( T ) = i aJ ) 2 COS {UjT + 'Ipi ) 2 
wo (ipl, , eine gleichverteilte Folge in E s ist. Weiter ist dann 


t J k= [ r k( ff ) da - 
Jo 


Wir erhalten die s Geschwindigkeiten 

v i = ^/tg^ + V’i) 
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und {ujki • • •) u sk) sind Lösungen des Gleichungssystems 


d 2 Uj 

dr 2 


T t üj Uj — 0 


mit den Anfangsbedingungen 

u jk (0) = a J cos (ip J k ), 


da 


( 8 ') 


dr 


= —c Oja J sin 




Wir haben bis jetzt den Fall betrachtet, daß alle Kj positiv sind. Wir 
untersuchen nun den Fall, daß alle Kj negativ sind (man könnte natürlich 
auch den Fall nehmen, daß nur ein Teil der Kj positiv sind, und die ande¬ 
ren Kj negativ sind. Dies wollen wir der Kürze halber nicht besprechen). 
Dann definieren wir die aj durch die Gleichung (Kj = Kj) 


K ! 


a j 


UJ 

J 


El 

K 


( 2 ") 


Wir nehmen jetzt die uj von der Gleichung 


d 2 Uj 

dr 2 


— ujJuj = 0 , 


also 


uj = cos(u;t + 'ipj) 

(wir benützen jetzt die hyperbolischen Funktionen cos, sin, tan). Es wird 
dann 

r J = (uj) 2 = (cos (ujt + ipj)) 2 

und 


Vj = uytanu;(r + ipj). 

Es wird also (77/ = uy (r + ipj)) 

V ] Kj ( 2 1 

— + — - uy f (tan rjj) + - - 

2 Vj \ (COS T]j) 


= UJj 


für 7 = 1,..., j. Hier haben wir benützt, daß cos 2 — sin 2 = 1 ist. 

Wir wollen die beiden Fälle: Fall I (alle Kj größer 0) und Fall II (alle 
Kj < 0) noch weiter untersuchen, indem wir auf die Fälle eine stereogra- 
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phische Abbildung ausüben. Es wird sich zeigen, daß wir im Falle I, in 
dem trigonometrische Funktionen bei der Darstellung der Lösungen 
verwendet wurden, jetzt hyperbolische Funktionen benützen werden 
müssen, dagegen ist es im Fall II gerade umgekehrt. 

Besprechen wir jetzt Fall I und kehren wir jetzt zum Fall (10) zurück. 
Es ist also 






7=1 


(ii) 


Wir setzen für j 




11 

Sil- 

1 

1 

VYv 

(12) 

Es ist Ä\ + 

+ Af = l.Wir setzen nun 



£ 27-1 + £ 2 / — y/EAje Sj 

(13,) 


ty- 1 - = y/EAje-*;, 

(13 2 ) 

also 

£ 2/-1 = y/EAjCos'ßj 

(14,) 

und 

£ 2 y = VE Aj sin 

(14z) 


Wir führen nun die Variablen x \,.. , X 2 S durch die (stereographische) 
Projektion in bezug auf das Hyperboloid (10) ein: 


jr^izj-A + £2 j Aje 0 ' 

**-' + ** = = ^ l-^sintf, 

VE 



(15,) 

(15a) 
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Setzen wir 


ßj = }J x 2j-\ ~ x 2j = 


A s sin ß s 


xy -1 = ßj cos ßj, x 2j — ßj sin ßj 


Es ist insbesondere 


cos ßj — 


_ Xy-\ 


Wir setzen nun 


Es ist dann 


Wir haben nun 


K = x y- 


-i = Ve 


A s cos ß s 
1 — A s sin ß s 


Es ist nun 


costf. 


/i y (l — ydj- sin t9j.) = \pEAj 


— VEA, = j. sin 


{ßj - VEAj) 2 = ßjA] sin ß s . 


cos 2 ß s — sin 2 ß s = 1, 


also erhalten wir für die rechte Seite von (20) nach (19) 


2 a 2 ( [ Z 




-1 )=Atf-ßA 2 ,. 


Setzen wir dies in (15) ein und quadrieren die rechte Seite von (20) aus, 
so erhalten wir 

+A]) - 2VeA jßj + (E - z 2 )Aj = 0. 


( 21 ) 
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Dividieren wir diese Gleichung durch den Koeffizienten 1 + A, von 
\Tj und ergänzen auf ein Quadrat, so erhalten wir 



AjVe 

1+Aj 


2 


A 2 ea j a ? 

1 +A] K (1 +Aj) 2 ' 


Setzen wir für \i, die Definition (16) ein, so erhalten wir für ] = 1,. ., s 


A ) 


1 + A? Z 



Ajy/E V 

1 + A 2 J 


EA]A> 

(1 +A 2 f 


( 22 ) 


Wir haben vor uns s Flächen (j = 1,..., s) im (X 2 /- 1 , X 2 /, ^)-Raum. 
Auf der Gesamtheit dieser Flächen, die nur in der ^-Koordinate 
Zusammenhängen, liegt die r-Kurve, die wir vorher studiert haben. Wir 
setzen 


6 / 


K, 


(aJcos(u)jT + i> J k )) 


2 ’ 


V 2 


Vj ufj tanculY + '/r'() 


V~2 


Die Flächen hängen noch durch die Gleichung 

A\+ +A] = \ (23) 

zusammen. 

Betrachten wir nun den Fall II 



Statt dem Hyperboloid in Fall I haben wir jetzt die Sphäre 

+C 2 = £- 

Wir setzen jetzt für j = 1,..., s 



Es gelten dann statt (l4i), (14 2 ) die Formeln 

£ 2,-1 = VEAjCosdj (140 


£ 2 / = VEAjsintfj. 


(14') 


und 
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Wir haben statt den hyperbolischen Funktionen cos und sin die 
trigonometrischen Funktionen. Die stereographische Projektion lautet 
jetzt (sie ist tatsächlich die klassische Projektion in bezug auf die Sphäre) 


/— A. cos r d; 

* 2,-1 = VE- J . 

J 1 — A s sin 

(15,') 

und 


/— A; sin 

x 2j = Ve- j - / . 

J 1 — A s sin ■ 

(1Ü) 

Die Definition (16) ist nun abzuändern 


+ 

7 

L^> 

II 

(16') 


(17), (18), (19) und (20) bleiben unverändert, es sind nur die hyperbo¬ 
lischen durch die trigonometrischen Funktionen zu ersetzen. Jetzt ist 
für die weitere Rechnung die Formel (*) durch die Formel 

cos 2 'dj + sin 2 , d s = 1 


zu ersetzen. Dies bedeutet, daß nun (21) so lautet 

ß 2 (\-A 2 )- 2^/EßjAj + A)(£ + E) = 0. (21') 

Wir setzen nun voraus, daß 

^ 2 < 1 (**) 

ist. Der Fall A 2 — 1 würde bedeuten, daß alle anderen A; gleich Null 
sind. 

Wir dividieren (21') durch 1 — A 2 , ergänzen wieder in /jy auf ein 
Quadrat und erhalten 

f n ajVev . aAI 

\ J 1~A?J 1-Af (1 +Af) 2 ' 


Setzen wir wieder für ßj die Formel (16') ein, so erhalten wir rTori 


( 



+ X 2j 


AjVe 

1 -A 2 


+ 


4 2 


1 -A‘ 


EA 2 A 2 s 

(1 


(22') 


Füry = l,. ., s im {xy- 1 , x 2 j, Q-Raum setzen wir 


^2/-l — 


Kj 


(Acos(üJjT + r ip J k )) 2 


, Q, = Q, tan(uyT + ?/Q). 
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Schaffen wir in (22) noch die Quadratwurzel weg, so erhalten wir 


: 2j-1 


A 2 

_ x 2 + _J A 
^ 1 +^ 2 ^ 


E.4 


1 

(i+^, 2 r 


4Evd : 


f 2 /-l 


— — 




Es sind dies also Flächen vierter Ordnung, so wie dies auch im Fall I, 
wo wir ja Tori vor uns haben, der Fall ist. Wir wollen aber trotzdem die 
Formeln hinschreiben: 



+ X 2j + 


A J 


1 


? 2 - EA‘ 


1 +Aj 

( 1 -Ä 2 ) 2 



+ x 2i) + 


4 EAj 


3 

Wir behandeln jetzt den zweidimensionalen Fall 

d 2 x Kx 
dt 2 r 2 

wo x der Vektor (x\ , ^ 2 ) u nd r — y.x 2 + x| ist. Die zugehörige Schwin¬ 
gungsgleichung ist 

S + “ 2 " = °- ( 2 ) 

Es sei a eine positive Zahl, dann gilt 

auj 2 = K. 


Es sei 

u — (u\,U 2 ) — (A cos(a jt + 'ip), Bsm(u)T + VO)- 
Wir setzen 


/— ^ T> /— • ^ 

= y^cos- , d = \Ja sm- • 

Es wird also Zl 2 + B 2 = 0 . Wir setzen noch 

E = CÜT + 'ip 


( 3 ) 


( 4 ) 


und nennen sie, wie in der Astronomie üblich, die exzentrische Ano¬ 
malie. Wir können also den Vektor u in der Form schreiben 


a E ot . E 

u = ya\ cos-cos—, sin —sin — 
2 2 ’ 22 


( 5 ) 
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Wir setzen im Anschluß an Levi-Civita 


2 2 

x^ = — ü / 2 


X2 = 1u\it2- ( 6 ) 

Um zu viele Indizes zu vermeiden, schreiben wir x\ = >r, X 2 =y und 
haben also (/ = \/—T) 

x + D' — (*i + ^ 2 ) 2 (7) 

Weiter setzen wir für die Zeit/ 

/ = J rdr. (8) 

r heißt natürlich wieder die fiktive Zeit. Es wird dann unter Benützung 
der Additionstheoreme der trigonometrischen Funktionen 

x = <?(cos a + cos E) (9) 

y = a sin a sin E (10) 


= \Jx 1 +j 2 = a( 1 + cos ol cos E). 


Für den Geschwindigkeitsvektor v erhalten wir nach (8) 

_ f dx dy\ / dx dy\dr (dx dy\ 

V dt ’ dt ) \ dr’ dr ) dt \ dr : dr ) / 


dx dy 

— = — au sin r, — — au) sin a cos E. 
dr dr 


Es wird also 


v = -(— sin E , sin a cos E) 

1 + cos a cos E 


und für das skalare Quadrat der Geschwindigkeit 


(1 + cos a cos E)' 


( — sin 2 E + sin 2 a cos 2 E). (13) 
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Wir betrachten auch gleich den Ausdruck 

dy dx 


\x, x\ 


x — = y- , 
dt dt 


(14) 


Man nennt diesen Ausdruck bekanntlich die Flächengeschwindigkeit. 
Wir erhalten nach kurzer Rechnung 

[x,x\ = 2au) sin ol. (15) 

Diese Flächengeschwindigkeit ist also konstant (zweites Keplersches 
Gesetz). Aus (9) und (10) ist sofort ersichtlich, daß hier eine Ellipse mit 
dem Parameter E vorliegt von der Gestalt 

2 / \ 2 


X 

-cos ol 

a 


+ 


y 


a sin ol 


1 . 


Es handelt sich also um eine Ellipse mit Mittelpunkt (^,0), wo 

e = a cos ol = V 'a 1 — b 2 (16) 

ist, die große Achse ist a, die Nebenachse b = a sin a. Der eine Brenn¬ 
punkt dieser Ellipse liegt im Nullpunkt (0,0) (die Sonne) (erstes 
Keplersches Gesetz). Es ist e die Exzentrizität und 


e = cos ol — - 
a 


(17) 


ist die numerische Exzentrizität. Um die bekannte Darstellung der 
Ellipse in Polarkoordinaten zu erhalten, müssen wir den Winkel $ (wahre 
Anomalie) durch 

cos ol + cos E 


cos $ = — = 

r 1 + cos ol cos E 


und 


y 

sin $ = - = 


sin a sin E 


r 1 + cos ol cos E 


einführen. Es wird 

1 — cos q cos $ 


d sin a 


1 + cos a cos E 


(18) 

(19) 

( 20 ) 


Man bezeichnet oft auch = tt — $ als die wahre Anomalie. Aus (20) 
folgt sofort die gesuchte Darstellung 


a sin ol 


a sin a 


r - 


P 


1 — COS Q COS $ 1 + COS Q COS 1 + COS OL COS <$' 
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Es ist 


p = a sin 2 ct = a{\ — e 2 ), 


der sogenannte Parameter der Ellipse. 

Es sei noch bemerkt, daß bei unserer Darstellung $ von der fiktiven 
Zeit r, ij> und co abhängt. Wenn wir jetzt die Flächengeschwindigkeit in 
Polarkoordinaten x = r cos $,j/ = r sin $ ausdrücken, so erhalten wir die 
berühmte Formel von Kepler 


t/<f> 

— = auo sin ct = laue 
dt 


* 


( 21 ) 


WO 

e* = Vl -e 2 (22) 

ist. 

Um in dieser Formelkette noch weiter fortzusetzen, die heutzutage 
trivial erscheint, die aber eine lange Entwicklung hinter sich hat, bilden 
wir uns 


und 


und erhalten 


1 + cos $ 


1 — cos $ = 


(1 + cos et) (1 + cos E) 
1 + cos a cos E 

(1 — cos et) (1 — cos E) 
1 + cos et cos E 


9 $ 1 — cos $ 1 — cos et 1 — cos E 1 — cos et 9 E 

tg — =-=-=-tg —, 

2 1 + cos $ 1 + cos et 1 + cos E 1 + cos et 2 

also die bis heute verwendete berühmte Formel 


$ /1 -e E 

tS 2 = y—e tg 2- 


(23) 


Wenn wir noch einmal auf (11) zurückkommen, so erhalten wir für 
E = 7t das Perihel Ii = a{\ — cos et), also wenn wir cos et = e positiv 
annehmen, den kleinsten Wert für r, also das Perihel, und für E = 0 den 
größten Wert, das Aphel. Es muß dann bei unserer Darstellung 

LOT + ijj = 0 

sein, also 



Wir werden darauf noch zu sprechen kommen. 
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Wir können aus (21) den Winkel $ noch explizit durch / darstellen: Es 
ist 


Es ist also 


also 


auj£* 

dt r 2 


d<& d$> dt au)£* 

dr dt dr r 


$=LÜE* 


dr 

1 + e cos E 


Wir wollen noch t als Funktion von E darstellen 


t = a (1 + cos <2 cos E)dcr, 

Jo 

wobei 


war. Es wird also 


E = ujg + 


t — a\ T 


£ sin E 


tu 


+ c, 


wo C eine Konstante ist, nämlich bei uns 

a sin 'ip 


C=- 


Cü 


Es besitzen cos E bzw. sin E die Periode (tt/u;). Es wird daher die 
Umlaufzeit 

a2ir 


T = 


LJ 


Nun ist cj = —, also wird 


,3/2 


T = ^2tt, 

Vk 


also 


T 2 _ 4tt 2 


(drittes Keplersches Gesetz). 
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Es ist nun (3) und (11) 


K 

r 


ÜJ 


1 + COS OL COS E 


(24) 


Benützen wir (18), indem wir noch die Formel sin 2 a = 1 — cos 2 ol, ber¬ 


ücksichtigen, so erhalten wir 

v 2 co 2 1 — cos 2 ol cos 2 E 


2 2 (1 + cos ol cos E) 2 

dann wird also 

v 2 K tu 2 (cos 2 ol cos 2 E — 1 + 2(1 + cos ol cos E) 
2 r (1 + cos ol cos E) 2 

also wieder die Energiegleichung 

ii 2 v 

.2 


(25) 


— = —üü 


(26) 


XX +jj = - —— 


(27) 

(28) 


Differenzieren wir (26) und die Gleichung (15) nach t, so erhalten wir 1 

Kdr_ 
r 2 dt 

yx — xj = 0. 

Die Gleichungsdeterminante ist nun, wenn r ^ 0 

dr 

xx + y y = r — 0 

dt 

und wir erhalten die Newtonschen Formeln 

.. K 

x = — x 

r 5 

.. K 

J = ( 29 ) 

Wir betrachten noch kurz den Fall der Hyperbel. Wir brauchen nur 
wieder die trigonometrischen durch die hyperbolischen Funktionen zu 
ersetzen. Wir erhalten 

x = ^(Cos ol + Cos E) (5 r ) 

y — a Sin a Sin E (6') 

r = a(\ + Cos ol Cos E). ( 7 f ) 


'Die Gleichung (28) ist auch an und für sich interessant und steht auch in anderer 
Weise in enger Verbindung mit der Schwingungsgleichung (vgl. P. Funk, Variationsrech¬ 
nung). 
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Es wird 




cos a 


J 


cos $ 


sin $ — 


a sin et 

cos OL + cos E x 
1 + cos a cos E r 

sincr + sinE _y 
1 + cos a cos E r 


Es wird 


1 — cos öl cos ip = — 


sin 2 ol 


Wir erhalten 


bzw. mit = 7t — <1? 


a( 1 + cos a cos E) 

P 


1 — COS OL COS $ 

p 

1 + COS OL COS $ ’ 


wo p——a sin a ist. 

Wir setzen noch e = cos a > 1 und e — ae. Es wird 


\x, x \ = a lü sin ol. 


(**) 

(18') 

(19') 

( 20 ') 

(+) 

(++) 

(15') 


Wir können nun, wie schon in §2, mehrere Planeten mit zugehörigen 
U\, ., co s mit 

cj 2 + + ujf = E 

betrachten (Pythagoras nannte dies eine Harmonie). 

Wir wollen noch studieren, was geschieht, wenn wir, sagen wir im 
elliptischen Fall, eine Drehung auf diese Kegelschnitte ausführen. Neh¬ 
men wir eine Drehung mit der Drehgeschwindigkeit H/. Es sei 

x = £ cos fit + ?7 sin fit 
j = — £ sin 0/ + rj cos fit. 

x — £ cos 0/ + f] sin fit + fly 
y — —£ sin fit + f] cos fit — fix. 


Es wird 
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Daraus folgt 


x 2 +j 2 = £ + f] 2 + Q 2 (x 2 +j/ 2 ) + 2Cl(xy —yx) 


und 


xy —yx = £t) — r)£ + D(x 2 +j 2 ). 
Wenden wir (18) und (26) an, so erhalten wir 

- + 7 ) 2 ) - Q 2 (£ 2 + Tj 2 ) - — ^ = -cü 2 - VLauj 


sin a 


und 


tv ~ = aüJ sin et — D(£ 2 + rj 2 ). 


Weiter ist jetzt 


r l — £ Z + ?? Z = ■ x ' 2 +J 2 


Wir differenzieren wieder und erhalten bei Auflösung der beiden 
Gleichungen 

s-2Mi=(^-n 2 y 

K 


jj + = ( — - Q ) r\. 


Führen wir auch in §1 eine Drehung durch 


so erhalten wir analog 


x = r cos 0/ 
y — r sin f2/, 


weiters 


x — r cos Qt — Qr sin 0/ 
y = r sin Qt + Dr cos Q/, 


X 2 +jj / 2 l.o 2\ K 
—-X = -(0 2 + fir 2 ) =- 
2 2 r 


CJ 


Dr' 


—yx + xy = Dr 


2 


und 
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Differenzieren dieser beiden Gleichungen und das Auflösen der differ¬ 
enzierten Gleichungen liefert 


2 2 2 
mit r = x +jy 


x - 2 Qy 

+ 2 fix 



X 

r 



J 

r 


Anhang zu §3 

Ersetzen wir in Formel (18) und (19) cos durch X', sin durch Y so 
erhalten wir nach G. Kowalewski eine Transformation, welche in der 
nichteuklidischen Geometrie wohl bekannt ist: 


X' 


X + cos a. 
cos aX + 1 ’ 


. 2Y sin a 

Y =-, 

cos olX + 1 


(*) 


welche den Kreis X 2 + Y 2 = 1 in den Kreis X /2 + Y ,Z — 1 überführt. 
Ersetzen wir in (*) die trigonometrischen Funktionen durch die hyper¬ 
bolischen Funktionen, so erhalten wir die Transformation 

. X + cos a . 2 Y sin a , . 

X =-, Y — -, (**) 

X cos a Y 1 X cos a + 1 


welche die Hyperbel X 2 —X 2 = lin den Kreis X' 2 — Y' 2 
Führen wir homogene Koordinaten 


1 überführt. 


X = -, Y = -„ bzw. X’ = * 


Y f = 


V 


C C ■ C" ‘ C 7 

ein, so erhalten wir im elliptischen Fall die projektive Transformation A. 

£ 7 = ^ + £ cos ot 
rj = 2r\ sin a 

C = £ cos ct + C- 

Ordnen wir der Transformation A. die Matrix A (et) 


A(ol) = 


1 

0 

COS OL 


0 COS OL 

sin ct 0 

0 1 
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mit Determinante D — sin 3 a zu, so erhalten wir 


A~' 


1 0 COS OL* 

0 sin a* 0 
sin ol* 0 1 


D _1 


mit ol* = — ol. 

Betrachten wir A (c^), A(ol 2 ), so erhalten wir für die zusammenge¬ 
setzte Transformation .^4 (et}), A(cl 2 ) die Transformation A(a 3 ) mit 

cos OL\ + cos cl 2 
cos ol 3 =- 

1 + COS 0L\ COS 0L 2 

sin ol\ sin ol 2 

sin ct 3 =-. 

1 + COS OL\ cos ol 2 

Wir haben also eine Gruppe vor uns. Setzen wir sin a = tan ß, so erhalten 
wir 

tan ß\ + tan ß 2 

tan P 3 =---— , 

1 + tan ß\ tan ß 2 

also 

ßi = ß\ + ßi- 

Im hyperbolischen Fall (**) ist Sin a durch tg ß zu ersetzen. 

Es liegt nahe, gleichverteilte Folgen (a N ), insbesondere solche, die zu 
pythagoräischen Tripel führen, und die zugehörigen Transformationen zu 
studieren, vgl. zu diesem Gegenstand auch §4. Wir wollen in einer ande¬ 
ren Arbeit darauf näher eingehen. 

Interessant ist die verkürzte Transformation 


r?i 


1 

COS CL 

Y 

,c\ 


COS OL 

1 

_C. 


bzw. 


b'l 


1 

cos a 

Y 

,C'_ 


COS OL 

1 

_C_ 


mit Determinante 1 — cos 2 a bzw. 1 — cos 2 a. Im zweiten (hyperbolischen) 
Fall wird durch die inverse Abbildung der Kreis £ /2 + ?7 /2 = 1 in 


£ — T] 2 = 1 transformiert. 
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4 

Wir führen nun auf §3 (9) und (10), welche wir jetzt in der Form schrei¬ 
ben, 


x = ö(cos et + cos E) 

(i) 

J = a sin et sin E 

(2) 

f*«i 

II 

o 

(3) 

eine Drehung D durch, welche wir als Produkt 


D = JD 1 D 2 D 3 

(4) 

schreiben, wo 


cos 0, sin 0 0 


D 1 = — sin Fi cos 0 

(5) 


[ 0 0 lj 




Gleichverteilung 


51 


Weiters besitzt das vektorielle Produkt 

[X, V] = [X,X] = {[YZ], [ZX\, [XY]) 

die Komponenten 

(cos fl sin i , sin fl sin i, cos i) [x, x ], (11) 

wo [5c, x] in Formel §3 (15) gegeben ist. Da Det (X, Y, Z) = 0 ist, so ist 
[YZ]X + [ZX] Y + [XY]Z = 0, 

also liegen X, Y, Z in der Ebene 

cos fl sin iX + sin £7 sin iY + cos iZ = 0. (12) 

Schneiden wir diese Ebene mit Z = 0, so erhalten wir für die Schnittge¬ 
rade (aufsteigender Knoten) 

f = -tg« (*) 

und, wie aus (11) folgt 

[YZ] 2 + [ZX] 2 + [XY] 2 = * 2 cu 2 sin 2 a. (13) 

Wir erhalten weiter die Gleichung (nach §3 (29)) 

d 2 X __ KX 
dt 2 r 3 

Wir können aus der Drehung D neue Drehungen 

= D’lD^D'l 

erhalten, wo jetzt statt fl, i, 7 kfl, li , nr y steht. Wählen wir fl, i, 7 als irra¬ 
tionale Zahlen, welche mod 1 linear unabhängig sind, so können wir, 
wenn wir die ganzen Zahlen k , /, m passend wählen, jede beliebige Dre¬ 
hung beliebig genau approximieren. Wählen wir noch dazu eine gleich¬ 
verteilte Folge ((/%,) bzw. (e k ), so können wir, wenn wir cu als fest und 
damit a als fest annehmen, jede Anfangslage des Systems beliebig approx¬ 
imieren. 

Es bleibt noch a festzulegen. Wir nehmen einen Punkt (x 0 , Jo, Jo) als 
Anfangsbedingung. Wir wählen die Koordinaten dieses Punktes in fol- 


(14) 
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gender Weise: Es seien a, 7 zwei Zahlen inj 2 (/), also zwei verschiedene 
komplexe rationale Zahlen a und 7 und es sei 

xq -\- ijo = a 2 — 7 2 und = 0:7 + cry, 


dann ist 


a = \J x o + jo + zo = M 2 + N 2 > 0 

und eine rationale Zahl. 

Wenn a oder 7 gleich Null ist, so ist — 0 und der Punkt liegt in der 
Ebene. 

Haben wir mehrere Planeten, sagen wir s Planeten, so wählen wir s 
Paare von Zahlen oy, 7 / -, wobei die 

a j = Wj\ 2 + | 7 y| 2 

verschieden sind. 

Damit sind auch die coj gegeben. Ob wir dadurch auch die Folge von 
Titus-Bode erhalten, möge bis auf weiteres dahingestellt werden. 

Natürlich können die Folgen auch zu einer beliebigen Dichte gleich¬ 
verteilt sein. 

Wenn wir nun auf die Dichte Eins spezialisieren, erscheint es angenehm 
z.B. solche Folgen zu nehmen, für die e = cos a, e* = sin a rationale Zah¬ 
len sind, d.h. pythagoräischeTripel vorliegen. Das gilt auch für 

E 

cos — 

2 

da ja 

A B 

cos (p = - - =, sin (p = , — 

\JA 2 + B 2 V A 2 + B 2 


= cos (ujt + !pj ) = cos <p, cos uJt + sin i pj sin uj/, 


war. 

Mit pythagoräischen Tripeln habe ich mich schon früher beschäftigt, 
so in meiner Arbeit ([ 1 ]) und in der Arbeit ([ 3 ]). 

2 Dieses Gesetz, welches vielleicht nur approximativ gültig ist, hat stets Interesse auf 
sich gezogen. Der Verfasser der vorliegenden Arbeit hat schon in seiner Jugend den 
Versuch gemacht, durch Interpolation der vorliegenden Daten, ein solches Gesetz zu 
finden. Es war wohl C.F. von Weizäcker der erste, der eine Herleitung dieses Gesetzes 
versucht hat (Z. f. Astrophysik 22 (1944)). Der Verfasser plant, die Gedankengänge von 
Weizäckers mit den vorliegenden Gedanken zu kombinieren. 
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Wir betrachten hier s Primzahlen p h ... ,p s von der Gestalt Ak +1, z.B. 
für s — 6 die Primzahlen 5,13,17, 29, 37,41 (Zahlenmystik). 

Jede solche Primzahl pj läßt sich im Gaußschen Zahlkörper — 1) 

eindeutig schreiben 

Pj ~ WW’ 

wo tt j und TTj Primzahlen inj2(\/— 1) sind. Wir definieren nun Zahlen v,- 
durch 


TT; 


TT; 


o'Xj 


für j — 1, ., s. Es sind diese s Zahlen dividiert durch 27t modulo 1 linear 

unabhängig. Es ist also die Folge (k = 1,...) 

&X\ &Xs 

2t r ’ ’ 2tt 


im Einheitswürfel E s gleichverteilt mit einer Diskrepanz (K absolute 
Konstante) 

K(logNY 


Es ist (L natürliche Zahl) 


D s < 

Un - N 1 A 


.V=iL Xj _ ( 

\tt j 


und wenn L gerade 2g ist und 

TTj — cij “b y /— 1 bj 

(< Tj , bj natürliche Zahlen) geschrieben wird, so wird 
ttJ = aj — bj + 2\J —\djbj 
\*j I = a i - b j ■ 


cos 2 gpj = 
sin 2 ,g<pj = 


2 - b 2 


if+ b I 

2 a.jbj 

7Tb 2 


a jg 


2 -b 2 


j& 


“}i + k M 


_ 2 a jgbjg 


a M + b M 


so ist 
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Wir wollen Q, /, 7 in der Formel für die Drehung D zusammen behan¬ 
deln. 

Die zugehörigen Primzahlen, welche die zugehörigen pythagor- 
äischen Tripel liefern, wollen wir mit p\,p 2 ,p?> bezeichnen, also z.B. 5,13, 
17. Dies gelingt so, daß wir die zugehörige unitäreTransformation 

*11 ayi 
*21 *22 


aufstellen, aus der dann in bekannterWeise die Darstellung von D folgt. 

Diese Darstellung von D durch pythagoräischeTripel gestattet weitere 
Anwendungen, die hier nicht mehr behandelt werden sollen. 

Wir setzen 


*11 


*21 


7T|7T2 <23 

MMM ’ 


7Ti7T2 ^3 

MMM' 


*12 — 


*22 


7Tl7T2 b 3 

MM M 

7Ti7T2 ^3 

MKIW 


(15) 


also z.B. wenn wir das Tripel 5,13,17 nehmen, so wäre 

77 — 1 T 2\f—\ 

7T 2 = 3 + 2 VM 
7T 3 = 1 + 4\/—T. 


Es ist 


<^ 11^22 — 0420*21 — 1 - 

Weiters setzen wir 

*0 = ~ (ot 11 + 0122) = (*1*2 — b \ ^2) *3 /* 2 
*1 = ~ (0*12 + 0(21) — (*1*2 + ^1 ^2)^3 /* 2 
*2 = — — («12 - «21) = (£1*2 - hMbi / a 2 

*3 — ~ “ (ü!ii — CÜ22) = (*1^2 + *2^1 )*3/* 2 ) 

WO 

ß 2 = ( ß 2 + ^ 2 )(* 2 2 + ^ 2 2 ) fe 2 + ^3 )• 


( 16 ) 
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Zu (1) gehört die unitäreTransformation 

also explizit 

Z' 


«li-Z + er 12 

(17) 

«21Z + Ct22 

7T2 l^\d 3 Z T i'Kjb 3 

7T2 ilt\b 3 Z T 71 "!<23 

(18) 


Zu (2) gehört das Quaternion 

{a\d2 — b \+ [d\d2 + b\b2)b 3 e\ 

+ {b\d2 ~ bid\)e2 + (<^ 2^1 + b\d2]d^e^, 

wo e 0) e h 6 2 , e 3 die bekannten Quaternioneneinheiten sind. 

Die Transformation (17) schreiben wir kurz 


« 11 

«12 


7Ti7T2^3 

i-KxTT 2 b 3 

.«21 

«22 _ 


_/7Ti7f 2 ^3 

Tt\Tt2dl 


Wir erhalten für die Drehung D die folgende Darstellung 


wobei 




C\\ 

C\2 

C 13 " 





D = 

Ql 

C 22 

6-23 

) 





_ Ql 

C 32 

C 33 _ 





Qi = 

2 1 

d 0 + 

2 
d | 

2 

d 3 d 

2 

3 




C 22 = 

2 

a 0 

2 1 

d^ + 

2 

d ly 

2 

3 




Q 3 — 

2 

a 0 

2 
d | 

2 1 

#2 T ö 

2 

3 



C 23 ~ 

- 2 {d2d 3 — 

dQd\) 

Q2 

— 2(^2^3 

+ 

dQd\) 

Ql - 

= 2(^3^! — 

dQd 2 ) 

C , 3 


+ 

dQd2) 

C\2 ~ 

= 2{d\d2 — 

dQd 3 ) 

Ql 

CM 

IS 

'S 

cT 

II 

+ 

dQd 3 ) 


(19) 


Man kann jetzt für (16) einsetzen. Eine ausführliche Darstellung findet 
sich in der Arbeit des Autors [4]. 

Man kann damit auch sphärische Dreiecke berechnen, dies soll an 
anderer Stelle ausgeführt werden. 
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Anhang zu §4 

Es liege eine Folge (x, ? ) modulo 1 vor, welche zur Dichte p mit Diskrepanz 
D P N gleichverteilt ist. Es sei L eine natürliche Zahl, dann betrachten wir 
nun die Gleichverteilung (wenn vorhanden) der Folge 


Bekanntlich gilt für jede Funktion F von beschränkter Variation V(F) 
im Einheitsintervall E 

— ^F(^)- / F(x)p{x)dx <V(F)D p n . 
iV k =l Jo 

Es sei nun e > 0 gegeben, welches wir erst später wählen werden. Es sei 
dann für £ < x < 1 

F M =f(-T~ [-Vh 


F M = 

für 0 < x < £. Dabei sei f periodisch mit der Periode 1 und beschränkter 
Variation. Wenden wir die obige Formel an, so haben wir zunächst das 
Integral 

J(e) = £ F(x)p(x)dx 

zu studieren. Es ist 

J(ß) = 

Wir setzen x~ L = £ als neue Variable, dann haben wir 

j(e)p (V ,/L ) r (,/L+1 V(«k 

=i?TX r,/i ) fM/<£V£ 

= zL' [ p(«+-*r ,/L )/(ü«+ü _(,/L+1 W 
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Setzen wir 


so ist also 


J(e)= ['f(S)Pi(SH- 

Jo 

Wir haben dabei einen Fehler gemacht, da ja k nicht bis ins Unendliche 
geht. Da aber p beschränkt ist, 0 < p <A und da (£ + k) 
konvergent ist, so ist der Fehler < £ AM, wo M(/) eine obere Schranke 
für /ist. 

Wir wollen nun noch V(f) abschätzen. Nehmen wir der Einfachheit 
halber an, daß f sogar stetig differenzierbar ist, dann ist im Intervall 

c < x < 1 

' W 

-L 




also wenn Mi eine obere Schranke von f' ist, 


,L+i 


f - L 

x 


< Mi 


rL+ i 


Wir haben also insgesamt 

iE Ai 


x„>£ 


- / f(x)p\{x)dx 

Jo 


< D p 
— u N 


Mi 

pL+l 


+ A + e(l + MA), 


da ja 


1 >r~-v 

N ? 1 

X„<£ 


p(x)d: 


x 


< eAD p n . 


Wählen wir nun = D^/^ L+2 ^, so erhalten wir 


— y^ f( — 

N J V x L 


e.v <x„ 


L 

L J 


xt 


— / f{x)p\{x)dx 

Jo 


< C(D» n ) 


P n(1/(L+2)) 


Der Fall L, = 1 wurde schon (ohne Restabschätzung) von Dirichlet 
behandelt (man vgl. meine Arbeit [3]), wo die Abschätzung auf anderem 
Weg gefunden wurde). 
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Analog behandeln wir den Fall 

1 


( sin | x /; ) 


1 


Hier setzt man 


(sinfx,,) 


( sin | x„) J 

r = e 


Die neue Dichte wird 


m 2 arcsin(£ +/ff) V ) d£ 

= / = 77T -v+ 1 / r- 

* \/i - (£+*r (2/L) (i+i:) 


5 


Es ist bemerkenswert, daß ein solcher Zusammenhang wie in §1 auch bei 
dem Friedmannschen Modell des Kosmos auftritt. Es handelt sich, wenn 
R der Weltradius ist, im elliptischen Fall um die Differentialgleichung 
(D>0) 


dR\ d 2 R 7 

Tr) = 2R ^ + WR ’ 


(1) 


wobei r die Eigenzeit ist, d.h. die Zeit, die durch eine vom Körper mit¬ 
geführte Uhr angezeigt wird, die in §1 fiktive Zeit genannt wurde. Die 
Zeit t wird wieder durch das Integral 

t = J R{r)dr (2) 


gegeben. Man nimmt üblicherweise u = 1. Wir denken uns die Eigenzeit 
durch cüt gegeben. Wir differenzieren nun (1) nach r (wir wollen jetzt die 
Ableitungen r durch einen. Punkt andeuten, z.B. statt ^ kurz R schrei¬ 
ben, dagegen wollen wir die Ableitungen nach der Zeit t in der üblichen 
Form R' statt ^ schreiben. 

Wir differenzieren nun (1) nach r und erhalten, da sich bei der Differe¬ 
ntiation die Glieder RR wegheben, 

R + w 2 R = 0. (3) 

Wenn wir R = v setzen, so erhalten wir für v die Differentialgleichung 

v + cj 2 u = 0. (3') 
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Die allgemeine Lösung ist 

v = A cos r + B sin r = C sin(urr + cp ). 
Integrieren wir (4), so erhalten wir, wenn wir 



Lü 


setzen, 


R = 


vdr = —K cos E + L, 


( 4 ) 

( 5 ) 

( 6 ) 


wo E = cur + cp ist. 

Aus (4), (5) und (1) folgt nun 

R = v = lj 2 K cos B = cu 2 R cos E. 

Setzen wir (1) ein, so erhalten wir nach (4), (5) und (7) 

K 2 lü 2 = (L - K cos E)u 2 (L - K cos E) + 2üj 2 K cos E 

= cj 2 (L — K cos E)(L + K cos E). 

Daraus folgt 

K 2 sin 2 E = L 2 — K 2 cos 2 E, 

also 

K 2 = E 2 , also E = ±K. 

Da R > 0 sein soll, so ist 

R = K(1 - cosE) = 2K sin 2 -. 

Es wird 



( 8 ) 


( 9 ) 






(10) 
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(10) nennt man die Friedmanngleichung. Wir differenzieren nun und 
erhalten 


also an allen Stellen, wo 
ist, erhalten wir 


/ // 2tu 2 , 

R R H-— R = 0, 

R z 

R' ± 0 


R" = - 


2uj 2 

R^’ 


also folgt aus (10) 


RR" = -?^=- w : 
R 


R' 


also folgt aus (9) und (10) 
RR" _ 
~R?~ ~ 


UJ^ 

BR 2 ~ 1 


1 


—UJ 


;+l 


«g 2 f 


UJ 


( 11 ) 


Für die Hubble-Konstante erhalten wir 

u;sin E 


sin 2 # 


E 

UJ COS^r 


R 

TT _ _ _____ 

ß2 r r(* t-\2 


(^2) 

E(1 — cos E) z 2Xsin 4 # 

Wir benötigen nun noch die Gleichung für die Dichte p nach der 
ersten Feldgleichung der allgemeinen Relativitätstheorie 

-KMV + (..-r 


UJ 


sin 


2 ) \ 2 

2 E (ctg 2 1 + sin 2 f ) \ _ uj 2 


sin 8 # 


sin 6 #' 


Es wird also die Dichte 


UJ 


sin 6 # 


p(E, (/?) 
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Es bleibt noch t zu berechnen 


f T . . ( sin(o;r + ip) \ 

Wir erhalten noch eine Keplersche Gleichung. 
Betrachten wir jetzt den hyperbolischen Fall 


R" = 2RR" - u 2 R 2 , 


so erhalten wir wieder mit V = R 


Es wird also 


v — = 0. 


v = ujK sin E 


R = K(cosE- 1) = 2K sin 2 — 


und die Hubble-Konstante 


RR" _ tu 2 
- sin 2 § 

ite 

u; sin E cj cos f 

K(cosE-l) 2 ~ 2K sin 4 f 


Für die Dichte erhalten wir 


tu" 

P ■ 6 F 
sin # 


Es wird die zugehörige Keplersche Gleichung 

t = k{ sin(cjr —d —-h 


Wir können wie in §1 Vorgehen und s Welten mit den Eigenzeiten 

UtiT, ,W S T 


+ + ojf — E 
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und Anfangsbedingungen gleich verteilter Folgen ..., (pf) mit 
Dichte p und Diskrepanz betrachten und erhalten so ein Modell für 
die Viel weltentheorie von D. W Sciama. 

Wir wollen zum Schluß noch auf die Keplergleichung, die wir schon 
mehrmals zu betrachten hatten, 

( sin(üJT + (f) 

t = a\ T - 

V ^ 

eingehen, die wir so schreiben können 

1 

- (c ot + <p>) = E — sin E, (17") 

a 

dabei haben wir, wie schon früher E = cur + p> gesetzt. 

Bekanntlich wird die Gleichung 

M = E — e sin E 



im bezug auf E als Funktion von M durch die Reihe 

OO 

E(M) -M = ^2 sin LM (*) 

fei L 


gelöst, wobei J L die ganzzahligen Besselfunktionen sind. Die Reihe ist 
für 0 < £ < 1 absolut konvergent. Für £ < 1 ist dies leicht einzusehen, der 
Fall e — 1 ist schwieriger. Es gilt zwar eine Aussage für JjJE), eine asymp¬ 
totische Formel, die ziemlich intrikat ist. Wir brauchen aber eine Abschät¬ 
zung nach oben. Dazu gehen wir von der Formel 

i r 

Jl(L) = — / cos E(a — sin a)da (18) 

TT Jo 


aus. Wir zerlegen das Integral in zwei Teile. Der erste Teil sei für ein 6 > 0 

rix— 6 

A— cos L(a — sin a)da 

Js 


und der zweite Teil sei 



Es ist unmittelbar klar, daß 


ist. 


|B| < 26 
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Wenn wir den Integrand in A in der Form schreiben 

(1 — cos a) cos L(a — sin ct) 

1 — cos a 

so liefert partielle Integration nach a sofort 
sin L(a — sin a) 


B = 


L( 1 — cos a ) 


* i r ~ 6 . . . . d ( 

H-/ sinlQ! — sin a )— — 

o EJs äa\\ 


COS OL 


da. 


Es ist nun zunächst der ausintegrierte Teil zu betrachten. Es ist 

? a ? 6 6 2 

1 — cos a = 2 sin — > 2 sin - > —, 

2 ~ 2 _ 2 

also ist er höchstens 


L6 2 


Das Integral ist 


1 

< — 
“ L 


’ 7 t — 6 


d 


1 


da V 1 — cos a 


da. 


Nun ist die Ableitung nach a im Integrationsintervall positiv, wir kön¬ 
nen also integrieren und erhalten die gleiche Abschätzung wie vorher. Es 
wird dann 


\JdL)\<46 + — . 


Wählen wir 6 3 = 1/L, dann ist 


L 3 


Es wird also 


I <20 ^ IA* < 40m' 


L—w 

Wir erhalten also 


L 


L—m 


1 Ä J l (L) 

E = cor + cp -) -- sin L{ujt + c p). 

L 


(19) 
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Brechen wir beim m-X&n Glied ab, so erhalten wir 

1 -J" > T (E) t? 

E = cot + 09 -——- sin L(tu/ + ip) H-f 

I, I f J T 1 


»bi L 


(M < 40). 


Wenden wir auf JjJE) eine Quadraturformel mit Hilfe einer 
teilten Folge (a k ) und Diskrepanz D* N an, so erhalten wir 

gleichver- 

Jl(L) = S W (L) + F, 

(20) 

WO 


| m 

S W (L) = — 7t cos (L(a k - sino^)) 

(21) 

und 


\F\ < LD‘ n 

(22) 


ist, da ja die Variation < jL. 

Setzen wir 

E:, = 2(w/ + V )--f' : ^7^smL(^ + v), (23) 

so ist 

E = El + F u (24) 

WO 

|Fi| < 40 + mD* n ). 


Wählen wir nun 


5 

o 

II 

1-1 

J-I^j 

1_1 

(25) 

so gilt für den Fehler F 2 


|f 2 | < 60 (D‘ n )\ 

(26) 

also haben wir 



E = E:, n + ^ 60 (Dy*. 
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Wollen wir z.B. r(/) aus §3 berechnen, so erhalten wir 

r{t ) = a(l ff- cos a cos E) = 1 + cos ct cos E * ff- F* 

mit 

F* < | cos E — cosF*i < |F-F*| < (D* N f 

Wir haben in §1 und [2] §7, indirekt aber in beiden Arbeiten, immer 
wieder Systeme Ei,..., £, betrachtet, wobei die Frequenzen co h ..cu s , 
welche eine Harmonie im Sinne von Pythagoras bzw. Kepler bilden. 
Diese Systeme werden durch Anfangsbedingungen, welche die Gestalt 
einer r-dimensionalen Folge haben, gleichverteilt bezüglich einer Dichte 
p, welche ein Produkt von Dichten p h ..., p s im bezug auf E 1? ., E x ist. 
Die Bestandteile des Systems Ei,..., Ej. stehen untereinander in keiner 
Wechselwirkung. Wir denken uns nämlich nach einer Idee von A. Ein¬ 
stein (in [8], S. 820 — 825), welche J. von Neumann in seinem Buch [5] 
abstrakt formuliert, jedes Ein einem Kasten (in einer Schachtel) Kj ein¬ 
geschlossen, welcher gegenüber äußeren Einwirkungen undurchlässig 
ist. Alle Fy werden in einem großen Kasten K eingeschlossen. Wir können 
diese E y als Moleküle eines Gases in K auffassen, so wie dies in der Physik 
bei der Mischung chemisch differenter Gase üblich ist, und dann die 
Gesetze der Thermodynamik anwenden. Übrigens können wir alle diese 
verschiedenen Systeme, welche wir in [2] §7 und §1 betrachtet haben, in 
einem Superkasten K zusammenfassen. 

Kehren wir zu unserem System E 1? . ., Ej. zurück und definieren wir 
als Energie des Systems 

E = tO\ ff- ff Lüf 

und als Entropie ff, welche wir als Mengenfunktion auf E! X x E ; 
auffassen, folgendermaßen: 

Sind A h ... ,A S quadrierbare Mengen in E t ,..., E J? so sei 
S(A U . ,A s ) = ^2 lg / Pj(xj)dxj. 

j= l 

Es ist also 

S(A U .. ,A S ) = log / / p 1 (x 1 )p s (x s )dx 1 ■dx s . 

JA\ Ja, 

Wir können K in ein Wärmebad mit Temperatur F geben, die freie 
Energie U — E — TS und dann die Gesetze der Thermodynamik anwen¬ 
den. 
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So sieht die Gleichverteilung des „Staubes“ im Friedmannmodell aus 

1 N 

cos $ 4 , cos A^, rjk(f) cos ^ sin A^, r& sin $^) — 


F(r cos <p cos A, r cos Lp sin A, r sin <&)r 2 dr sin $ \d<&d\ 


< V(F)D 6 n . 

F ist von beschränkter Variation. 


Anmerkungen 
Anmerkung 1 zu §1 Seite 14 

Betrachten wir noch den Fall tu = 0, also die Gleichung 


mit der allgemeinen Lösung 


= a(r + 'ip) 


und setzen 


dann wird 


Es wird 


x = u 2 — 1 = (a(r + 'ip)) 2 — 1 
y = 2u — 2 a(r + ip), 


r = u 2 + 1 = (ä(t + ip)) 2 + 1 . 


t = / rVr = 


! (r+ -0) : 


+ r + ip + C 


Wir haben also eine Parabel 

j 2 = 4(x+ 1). 

Es wird x — 2a 2 [r + ip ), y = 2a, v = - 


Es wird 


wenn wir K — 2a nehmen. 


v 2 = 4 ^ 2 


v 2 K 

-= 0 , 

2 r 
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Wir setzen 


u 2 — 1 . 2 u 

cos $ = —-r, sin <P = 


u 2 + 1 


1 + 


1 + cos $ 


1 — cos $ 


2 


also 


ü/ 2 + 1 

2 

* 2 + r 

2 


r — 


1 — cos $ 


Wir können die Keplergleichung explizit lösen. Wir haben die 
Gleichung 

w 2 ’ + — 3 t = 0 

zu lösen. Die Diskriminante dieser Gleichung ist 

—4 • 3 3 — 27 1 2 < 0, 

also nach Cardano 


w = 


\ 2 


/ 

- + 


+ 1 + 


\ 


~~2 + \{~2 1 + 1 ' 


Setzen wir 


~ 2 + \{~2 1 + 1= ct g “ 


+ 1 = — ctga, 


so wird 


w = - ctg 2 a. 


Anmerkung 2 

Eine Bemerkung zur Energiegleichung, siehe z.B. §1 Formel (6), 

v 2 2c ü 2 2 

= -2w 2 , 


2 cos 2 E 
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die eng mit der Schwingungsgleichung ü + cü 2 u = 0 mit E = lüt + p 
zusammenhängt. Wir können v = tü 2 als Frequenz und A = ctg(wr + ip) 
im Sinne von Broglie als Wellenlänge bezeichnen. 

Setzen wir S = log cos E + C, wo C eine additive Konstante ist. Wir 
können dies auch schreiben 


so ist 


S = log^d cos -E, 


dS sin E 

dr cos E 


und die Energiegleichung schreibt sich 


1 fdS\ 2 

2 V<9r ) 


K 

- = o, 

r 


wo K = 2üü 2 a und r — a cos 2 E ist. 

Führen wir eine zusätzliche Variable a ein und setzen 


W = — cu 2 a + S(t) = — vol + A(t), 

so erhalten wir 

dlV fdlE\ 2 K_ q 

da \ dr ) r 

Es ist 



dr 


im Sinne der Hamilton-Jacobischen Gleichung. 

Setzen wir mit Schrödinger, der diesen Weg später nicht mehr weiter¬ 
verfolgt hat, 

ip = K cos E, 

so wird 

dS i[) 
dr iß ’ 

dann wird die Energiegleichung 

Tp 2 — 2u 2 tg 2 E'ip z = 0. 
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Betrachten wir die Funktion F = ip 2 — 2cu 2 tg 2 Eip 2 als Funktion ip, ip, r 
und stellen die zugehörige Eulersche Gleichung auf, so erhalten wir für 
die Extremale Ydie Gleichung 

FI- + 2« Ve</<(t) = 0, 

für die man noch Randbedingungen vorschreiben kann. Wir erhalten 
eine Hillsche Gleichung. 


Anmerkung 3 zu §3 und §4 

Die Drehungen der Erde um ihre Achse sind ein Beispiel für eine Statis¬ 
tik von Drehungen. Nehmen wir in nullter Näherung für die Erde an, sie 
wäre ein kräftefreier Drehkreisel, dann können wir die Eulerschen Glei¬ 
chungen (A, B, C dieTrägheitsmomente der Erde) aufschreiben 

A d 4+(B- C)rq = 0, B^ + (C — A)rp = 0 
dt dt 

und 

^~dt + ^ ~ ^4 = 

Nehmen wir weiter an, daß A =B ist (auch dies ist näherungsweise 
erfüllt), so folgt aus der letzten Gleichung, daß r(t) eine Konstante, sagen 
wir r 0 , ist. Setzen wir 

A-C 
»=~— 

und 


27rA = jiro, 

so erhalten wir aus den beiden ersten Gleichungen 

— + 2ir\q = 0, — - 2ttA p = 0. 

dt dt 

Wir kommen also wieder auf die Schwingungsgleichung zurück und 
erhalten als Lösung 

p(t, = a sin 27t(A t + </?), 
q{t , a,(ß) = a cos 2ir(\t + Lp ), 


wo a und Lp Integrationskonstante sind. 
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Es ist 

.2 , 2 _ 2 

p + ^ — a 

Wir nehmen a als positiv an. Es ist 

Ü2 = p + ^ + r — a +r 0 
die Energiegleichung. 

Wir nehmen nun eine Folge (a k , <p k ) im Einheitsquadrat E 2 gleichver¬ 
teilt mit der Diskrepanz an. 

Es sei g eine Funktion von beschränkter Variation V(g) in E 2 und wir 
setzen 

G{a,(p) =g(a cos 27t(p, a sin 27np). 

Es sei 


N 


^N{g ) = — ^ G(a k cos 2 7T(p k , a k sin 27 up). 


k=\ 


Es ist 


wo 


An (G) = J + tfV{G)D N , 


J = G{a , ip)dadp> 

Je 2 

(|$| < l).Wir setzen a = a 2 (a > 0), also ist 

J = 2 G(a 2 ,Lp)oidad(p. 

Je 2 

Setzen wir 

£ — a cos 27np, r] — a sin 27rep, 

so wird 

/ = “ ff 

Führen wir die Drehung 

£ = £' cos 27tA t + rj sin InXt 
rj = —rj sin 27 tA/ + rj cos 2tt\ t 


(*) 
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im Integral J durch, so erhalten wir 

J 


£ /2 + t / 2 <1 


Setzen wir 


JJ cos 27 tA t + rj sin 27 tA/, 

— ^ sin 2tt\ t + rj cos 27tA t)d£jdrj 

= a sin 2mp, ?/= G!sin27r<£>, 


so erhalten wir 

J= ä u cos 27r( A/ + p\a sin 2ir(\t + p)dadp. 

Je 2 

Wenn wir nun von J zu 

1 N 

G(j>(t,a k ,<p k ),q(t,a k ,(p k ))=—'^2g(a k cos27r(\t-<p k ) : asm27r(\t + <p k )) 


k =1 


zurückgehen, so erhalten wir 


1 

N 


N 


^ ^ G[p(t 1 fi-ki Pk) ) tfij 1 ^ki Pk )) 


k=\ 


N 


^ ^ G(p(0i fi-ki Pk) 1 ^(0) ßki Pk)) 


< 2tiV(G)D N . 


So sieht bei uns der Liouvillesche Satz aus. 

Nehmen wir an, daß g die Indikatorfunktion i einer quadrierbaren 
Menge A mit der Schwankung cr(A , D^/ 3 '*) ist, so ändert sich nichts an 
der Überlegung. Es bleibt die Häufigkeit 

1 N 

Ty ^ ^ t(Aipi^tPk), q{t 5 a kiPk) s ) 

^ k=\ 


gleich der Häufigkeit dieser Punkte p(t, au, p>P) zur Zeit / = 0 bis auf einen 
Fehler von der Größe a(A, D ^ /3) ). 


Die Umlaufszeit 


27t 




ßr 0 


bleibt ungeändert. 
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Wir können die gleiche Methode auch auf die Planetenmodelle von 
Levi-Civita anwenden. Wir nehmen jetzt unter gleichen Voraussetzungen 

für g 

G(a , ip) — g{a sin 2 27 vcp, coa sin 27 Tip) 

(oj > 0). Es sei wieder 

1 N 


k=\ 


Wir führen wieder die Drehung (*) durch, allerdings schreiben wir r 
statt t und wir erhalten (u) — 2 ttX) 


J = 


Jf g((a cos 2 tt(\t + </?)) 2 , oja sin 27t(Ar + < p))dad(p. 


Jetzt kommt die Änderung: Wir setzen 

r(r) = <7(cos(u;t + 2tti p)) 2 

und 

r(r)v(r) = u;sin(a;r + 2tt<p). 

Es ist (r(r) < a < 1) 

v(t) = oJtg(u;r + 27 x<p). 

Wir setzen 

r k = a k cos 27 r<p k , v k = u) tg (W + 27 up k ). 

Es ist 

r 2 M + d) = 4<4- 

Wir erhalten also 


J 


= J'J^g{r(r,a,(p),r(T,a,ip),v(T,a, (p)dad<p 


und der Liouvillesche Satz wird 


N 


1 X JV A 1 N 

Jj^g( r ( T ,ak,Vk),r(T,a k ,(p k )v(T,a k ,(p k )^ g{r k ,r k v k ) 
k =l JN/ A=i 

<2K(G)D n . 


Dies bleibt natürlich richtig, wenn G nur quadratisch integrierbar ist. 
Es ist dann die Variation V(G) durch die Schwankung cr(g, D^/ 3 ^) zu 
ersetzen. 
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Ist i(K, R) die Indikatorfunktion eines Kreises vom Radius R in E 2 
mit Mittelpunkt 0, so ist 

1 N 

TT ^ ^ ['(K., t fc, tkVfe) 

k=\ 


die Häufigkeit der Punkte (qr k , r k v&), die zur fiktiven Zeit / = v im Kreise K 
liegen und dann ist sie auch zur Zeit r(/) bis auf einen Faktor, der von a 
und Djss abhängt, die gleiche. 

Wir haben stets in iE 2 : 0 < £ < 1, 0 < 77 < 1 gearbeitet. Alles geht 
analog, wenn wir ein anderes, z.B. ein größeres Rechteck, nehmen und 
wir können dies sofort auf E 2 zurückführen. 

Nehmen wir an, es liegen mehrere Kreisel oder Planeten oder mehrere 
Oszillatoren vor, sagen wir s solcher Objekte, so haben wir s Frequenzen 
Ai, ., Aj.. Wir haben Ns Anfangsbedingungen, die eine Zr-dimensio- 
nale Mannigfaltigkeit bilden. Wir wollen sie z.B. als gleichverteilte Folge 
in E 2t mit einer Diskrepanz D 2 ( ansehen, dann können wir alles übertra¬ 
gen- . 

Wir führen die Integration von Gleichungen des kräftefreien symme¬ 
trischen Kreisels zu Ende. Dazu müssen wir die Eulerschen Winkel er, ip, 
d einführen (der Eulersche Winkel er wird traditionell mit Lp bezeichnet. 
Wir haben aber bereits bei der Integration der Gleichungen für die 
Impulse p und q als Integrationskonstante mit Lp bezeichnet). Die Glei¬ 
chungen für die Winkel lauten bekanntlich: 


dd n dp) 

— = p cos a — q sin a, sin v —— = p sin a — q cos o 
dt dt 


da 

dt 


d'ip 

dt 


cos 




Wir tragen nun die bereits berechneten p, q, r ein: 

p — acos(Xt + (p), q — sin(A/ + Lp), r = tq, 

wo A = (1/27 t)(1 — CjA) ist (pp, i'o Integrationskonstanten). Zunächst ist 
= 0 Lösung. Im Falle 0 < d < (7 t/2 ) erhalten wir 

dd d'ip ( C \ da C 

dt ’ dt V A dt A ° 


also als Lösungen mit den Integrationskonstanten cr 0 , $ 0 , ^0 
$ = flo, cr = £/ + cr 0 , lp = 6t + 1p 0 . 
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Dabei ist 6 = (C/^4)(r 0 /cos f? 0 ) die Winkelgeschwindigkeit, mit der die Fi¬ 
gurenachse um die Impulsachse rotiert, also die Präzisionsgeschwindig¬ 
keit und 

£ = 1 + ( 1 -f) ro - 


£ ist die Eigendrehung, also die Drehgeschwindigkeit um die Figuren¬ 
achse. 

Die Drehung D wird gegeben durch 



D — D\D 2 D3, 



cos er, 

— sin er, 

0 ' 

D i = 

sin <7, 

COS (7, 

0 


0, 

o, 

1 


cos$, 

0, — sin 

d 

d 2 = 

o, 

1 , 0 



sin 7?, 

0, cos 

7? 


COS "0, 

— sin i / j , 

0 

d 3 = 

sin t/>, 

cos t/>, 

0 


0, 

o, 

1 


ist. Wir betreiben nun Statistik im Einheitswürfel E 6 mit Hilfe einer 
sechsdimensionalen gleichverteilten Folge 


( 2 

l dkl 

V * 


7? 


v 0/ Vf '/: \ 

“’27’ 2tt F 


Man wird, wie bei der Planetenbegwegung, die Winkel so wählen, daß 
wir pythagoräische Tripel erhalten. Besonders interessant ist der Kugel¬ 
kreisel A = B=C. Hier wird A = e = 0. Ein Beispiel ist der Würfel oder 
die Kugel. 

Wir können mit den entwickelten Formeln auch den Fall mehrerer 
voneinander unabhängiger Kreisel, z.B. Wü rfel, die nur durch Anfangs¬ 
bedingungen miteinander verknüpft sind, behandeln, so wie wir dies 
schon bei den Planeten getan haben. 


Anmerkung 4 zu §5 

Wir können nach Friedmann-Lemaitre-Milne-Heckmann mit den For¬ 
meln ein Weltmodell aufbauen. 
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Betrachten wir (die Striche sind die Ableitungen nach der Zeit) 

R' 

f ~ *R ’ 

dann folgt durch Differentiation 

. R" {R'\ 2 

f =Y- R 


also ist 


Wir setzen 


J ^J R3 


und es wird also 


pR 3 = K. 


Wir nennen p Dichte im physikalischen Sinn. Es wird 

^ = -3^R' = -3p-, 
dt R 4 R 


Wir haben also die Gleichung 

dp 


+ 3 pf — 0. 


Wir definieren nun Teilchen x — (x,j, z) 

x(t) = x 0 r(t), j(/) =7 0 r(/), ^(/) = zor(f) 
und setzen x = xor(t), (x 0 ,j 0 , zo) sei ein fester Punkt in R 3 Es sei 
_ (dx dy d%\ , x 

v ={*•*’*) = {Vu '*’ V3) 
die Geschwindigkeit dieser Teilchen. Es ist 

v = xf 
div v = 3/ 


Es wird also 


dp _ 

tt- + pdiv v — 0. 
dt 
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Wir haben also eine Kontinuitätsgleichung. Wir haben nun 

dpi) f dv\ d(viv 2 ) d^p A __ 
dt P \dxi dx 2 dx?, ) 

dp dvi ( dv 1 (dv 2 dv 3 

= Vi di + p ~di + p { 2v 'd7d v '\d^ + d^ 

= (-3/ 2 +/' + 4/ 2 Ei = -2px,. 

Es gelten analoge Formeln auch für j— 2 bzw. j — 3. 

Setzen wir 


$ — ^x-j + x 2 + x 3 J p, 

so haben wir (für j — 1, 2, 3) 

dp/) / dpQ dp 2 ) dfep h = d$ 

dt P \ dx 1 <9x 2 dxß / dx y 

Es sind also die x y (t) die Bahnkurven einer Strömung, auf die die Kraft 
mit dem Potential $ wirkt. Dieses Potential erfüllt die Poissongleichung 

= 6p. 

Dies rechtfertigt den Namen Dichte für p. 

Es ist $ bis auf Konstante das Potential einer Kugel K mit Mittelpunkt 
0 und Radius R und konstanter Dichte p in einem Punkt (x, J4 Z) im 
Innern der Kugel gleich 


V R ~ 3 ^ +J +Z ^ P )' 

Wir wollen nun zur Dichte p eine Folge (£^, p k ) konstruieren. Zunächst 
wollen wir sie normieren. Es sei 


p(r, <p)d<p = tp) 


und es sei 


Ppp = 


P(t, <p) 

VP) 


Wir bemerken noch, daß für r > 0 


P, V>) < 


ist. 
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Es sei (x k ,j k ) eine gleichverteilte Folge mit Diskrepanz D N , dann sei 


N r 


1 IN 

/u = -Y 

1 N f 


1 + X/ - / p*(ip)dip 

Jo 


1 +y k - 3 


P'/ 


rdr 


x /=l l 

Es ist die Diskrepanz dieser Folge 


< 


1 


(W)' 


■/(r) D 


1 n l/5 


N 


Es ist also für jede Funktion/ von beschränkter Variation V(f) 

= 3 / r 2 dr f(p,r)p*(T,ip)dT + Fehler, 
iV A=l 70 70 


WO 

IFehler| < K(/) ((u;r) 6 /(r)) _1 D^ 5 
ist. 

Es ist z.B. die Häufigkeit der Punkte (<£/, p k \ wo ^ < A, Tj k < R(r) 
(also die Häufigkeit des Staubes in der Kugel vom Radius R(r) 



Der Faktor 47t kommt von der Oberfläche der Kugel. 


Anmerkung 5 zu §5 

In der allgemeinen Relativitätstheorie kommt man auf das Fried- 
mannsche Modell in folgender Weise: 

Wir betrachten die Riemannsche Metrik 

ds 2 = —r 2 (t)dcr 2 + dt 2 , 


3 

da 2 = ^ ^ 7 jjdxjdxj 
v=i 


wo 
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ist und der Energietensor Ty mit den Komponenten 

Ty = 0 für 1 < i < 3, 1 <j < 3 

und T /4 = 0 für 1 < / < 3 und T 44 = p. 

Sollen die Einsteinschen Feldgleichungen erfüllt werden, so muß der ver¬ 
jüngte Krümmungstensor von der Gestalt 

Gy — 2 £ 

sein mit e = 0 , =h 1 und wir erhalten 

r' 

p + 3 — p = 0 

r 

und 

2e — 2r — rr — —r p. 
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